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Abstract 

We use the high statistics E-772 data on the nuclear dependence of the pro- 
duction of quarkonia {J ftp and T) and dimuons at large transverse momentum 
{pt) in p-A collisions to get information about the gluonic EMC effect. We 
find a satisfactory quantitative agreement of the theoretical predictions with 
the data although none of the models of the EMC effect we consider could 
account for the entire data. Since all the qualitative features are understood 
none the less in terms of perturbative QCD with nuclear dependent parton 
densities, our results suggest that these data can now be used for a better 
determination of the nuclear parton densities. Our conslusions are shown to 

be insensitive to the hadronisation mechanism for the quarkonia. 
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I. INTRODUCTION 



The observation of a nontrivial nuclear dependence of the parton densities, the EMC 
effect 0, confirmed later in a further series of deep inelastic scattering experiments 0, still 
lacks a clear theoretical understanding even after a decade. A large variety of models 0, 
with widely differing basic underlying physics mechanisms, have been proposed to explain 
the EMC effect. In all of these, some model parameters have to be fitted to reproduce the 
measured ratio p of the structure function Fg"^ (per nucleon) for the nucleus to that of the 
nucleon ^2^. As a result, the nuclear quark densities are very similar in all these models. 
However, their predictions for nuclear gluon densities are quite different. Hence a good 
measurement of the gluonic ratio. 



Pg 



AgP{x,Q 
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can play an important role in distinguishing the different models of EMC effect. This 
information can be obtained by studying the nuclear dependence of different hard scattering 
processes ||^ such as jet production direct photon production ||^, quarkonium {J/ip 
and T) production P,[7|J^,pil[] , ji'^jj,^ pair production [|TT|] with nuclear targets. Recently the 



associated production of J /if) and a photon at large pr has also been suggested as a 
means to catch a glimpse of the gluonic EMC effect and to help unravel the correct model 
of EMC effect. 

Nuclear dependence of quarkonium and p^ pT pair production can be very crucial for the 
various experimental signals of the formation of Quark-Gluon Plasma (QGP) in relativistic 
heavy ion collisions. It has been suggested [1^ that due to formation of QGP, the J/ip 



production in A- A collisions will be suppressed and the suppression is expected to exhibit 
a characteristic pr-dependence with a maximum suppression at p-p = |[T^. To analyse the 
feasibility of using this suppression as a signal for QGP formation it is, however, necessary 
to understand clearly the expectations of perturbative quantum chromodynamics, (pQCD), 
for the quarkonia production in A- A collisions. The observation of a pr-dependent J/ip- 



suppression by the NA38 experiment [|T5l further made such a study of pQCD imperative. 



The applicability of pQCD for this process is justified because of the short time scale for 
the production of a heavy quark-antiquark {QQ) pair which is 0{1/2Mq) for a heavy quark 
of mass Mq. In the pQCD approach the initial state dependence of J/ip ot T production 
comes through the nuclear parton densities. Earlier investigations have revealed that 

the nuclear dependence of parton densities gives rise to a J/ -i/^-suppression in A- A and p-A 
collisions which is qualitatively similar to the proposed QGP signature. Since any suppres- 
sion due to QGP formation or a hot dense hadron gas amounts merely to a change in the 
hadronisation model and/or final state interactions, it is clear that a critical assessment of 
the utility of J/'j/'-suppression as a signal for QGP formation needs a good understanding 
of nuclear parton densities. Since quarkonia are usually detected through their decay, 
the continuum pair production in A-A collisions has to be investigated as well. For- 

tunately, high statistics data on Jjtp oi T and pair production with nuclear targets 

ITtHIQI has recently become available from the Fermilab E772 experiment. In this paper we 
compute large-pr J/'^,'^ and dimuon production in p-A collisions in pQCD and compare 
it with the E772 data with a view to get a handle on the gluonic EMC effect. The main 
advantage of comparing the model predictions with all these data from the same experiment 
is of course in minimizing the effects of systematic uncertainties. We consider the processes 

p + A^ J/^(T,/xV) +^ 

as due to a hard collision between partons from the proton and the nuclear target with mass 
number A. The nuclear dependence is incorporated only through the parton densities which 
is where the different models of EMC effect come into picture. The hard scattering cross- 
sections for QQ production {Q = c (b) for J/iIj{T)) and fi'^fi~ pair production are available 



in literature p0|^2|. Two different hadronisation models p3| , p^ exist in the literature for 
the hadronisation of the QQ pair into a quarkonium. In our work we use both of them to 
assess the uncertainty in our results due to these. Earlier an attempt was made to extract 
the nuclear glue from the E772 data by using the approximation of gluon dominance of the 
partons in the initial state and by assuming that the QQ production takes place at pt = 



for the pair. We comment on the rehabihty of this approach. 

The plan of the paper is as follows: in the next section we summarize different theoretical 
models for nuclear parton densities (i.e. the EMC effect) which we use. In section III, we first 
present the details of the theoretical formalism we use to calculate the quarkonia production 
at large-pr. We then discuss and compare the two hadronisation models we employ and point 
out the major kinematical differences in the two models. Following this the limitations of the 
gluon dominance approximation, used in Ref. |^ to extract pg from J/ip and T production 
are highlighted. Finally, we end the section by presenting a comparison of our pQCD 
calculations with the data from the E772 experiment for the J / ip and T-production at large 
Pt- The theoretical predictions are obtained for the different nuclear parton densities by 
imposing the experimental cuts of the E772 group on the ranges of the kinematical variables. 
We also compare the predictions of different hadronisation models with each other. In the 
next section, we discuss the formalism for calculating the background, continuum /i^/i" 
pair production at large-p^ and present our results for the ratio of partially integrated cross- 
sections of dimuon production for different nuclear targets and parton densities and compare 
them with the E772 data. We then present our conclusions in the last section. 

II. MODELS FOR NUCLEAR PARTON DENSITIES 

We define the nuclear parton densities by removing a factor A, where A is the atomic 
number of the nucleus, so that the cross section defined below is per nucleon of the target, 
making it directly comparable to the pp case. The presence of the EMC-effect tells us that 
the nuclear parton densities are different from those of a nucleon even after this rescaling, 
i.e., the ratio p, defined in the introduction, is a nontrivial function of x, where x is fraction 
of the momentum of a nucleon which the parton carries: < x < A for the nucleus. There 
are many models in the literature which predict p in terms of some parameters which can be 
fixed by comparing with the data on the EMC-effect. The corresponding pg is then usually 
predicted in these models. We limit the scope of our present investigations by concentrating 
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on three different models, which are chosen as typical examples of sets of models based on 



similar physical concepts: i) the gas model ii) the rescaling model ||2^ and iii) the 
six-quark cluster model In all the three cases, we neglect Q CD-evolution, since such 
corrections are small in the kinematic range of interest to us. Fig. |1| displays the predicted 
Pg for all these three models. The details of these predictions are given below for each of 
the model. 



A. The Gas Model 



The parton densities of a nucleus with atomic number A are defined in the Gas model 



1 2^ as a sum of two components. 



^ r=l 



(1) 



The first term, occuring with a weight (1 — tu), is that for a free nucleon parton density 
after corrections for its Fermi motion inside the nucleus. The second component is written 
in terms of thermal distributions of momenta at a temperature T, leading to the following 
functions ff^{x] /i, T): 

1 



TT 



£n(l + ze-^) + -In^il + ze'^) + Ligl 

2 \z + el- 



(2) 



and 



TT 



(3) 



Here tq = 1.2 fm, = Mx/2T, z = exp(— /i/T), where M is the nucleon mass and Li2{x) 
is the Euler dilogarithm function. Using the constraint on the total baryon number of the 
nucleus to eliminate the chemical potential /x, one has two model parameters, T and u, for 
each nucleus. 

Using the CDHS parametrisations [P^], 
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F^P{x) = 1.1(1 + 3.7x)(l-x)3-^ 

x(Tp{x) = 0.17(1 -x)^-^^ (4) 
xfg/p{x) = 2.62(1 + 3.5x)(l -x)^-^ 

for -F^x) (here <7p{x) is the total sea density) and the data on p(x) = F^{x)/AF2{x) 
these parameters have been fixed [^] for many nuclei, including the ones used for the E772 
experiment. The corresponding Pg{x) is then predicted uniquely. We summarise them in 
Table |I[ It may be mentioned here that using a different set of structure functions for proton 
instead of eq. (^ necessitates a time-consuming re-analysis of the EMC-data to obtain T and 
u. For this reason we have not used more recent parametrisations for the proton structure 
function and we are also constrained to use different proton structure functions in different 
models. 



B. The Rescaling Model 



Scaling models p6|j29|j30[| seek to explain the EMC effect as arising from the change in 
the QCD-scale in going from a nucleon to a nucleus. The nuclear parton densities at a scale 

are obtained from the parton densities in a proton at the same by evolving them to 
a scale ^aQ^, i-e., the nuclear parton density per nucleon fi/AijiQ"^) is given by 



fi/A{x,Q'^) = fi/p{x,^AQ^) 



(5) 



In this paper, we use the rescaled nuclear densities as obtained in Refs. [ P^pOl where ^a = 
J^l"^ and the starting nucleon parton densitites were taken [^] to be a parameterisation of 
the EMC Deuterium data at = 20 GeV^. For further details, we refer the reader to Ref. 
0. 
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C. The six-quark cluster model 



The six-quark cluster model |27| is also a representative of the two-component models 
for EMC effect. In this model, it is assumed that when two nucleons get closer to each other 
than a certain critical radius they merge together to form a six-quark cluster. By assuming 
the probability to form higher clusters to be negligible, the remaining model inputs are 
the probability of forming such a cluster and the form of the parton distributions in the 
3-quark and 6-quark clusters. The latter are chosen using the quark-counting rules and the 
constraints of i) normalization of valence densities (an N-quark cluster has N valence quarks), 
and ii) conservation of momentum. It is further assumed that the average momentum carried 
by the sea partons is the same for the three- and six-quark clusters and that it is ~ 0.2 of 
that of the gluons. The forms of nuclear densities per nucleon in this model are, 

/,/^(x) = (l-e)/,3(x)+ I As (I) , (6) 



where e is the probability to find a six quark cluster which increases with A [31| and the 



subscripts denote the cluster size. The values of e which we used are given in Table |. 



Specific choices [|2^ for the valence density V{x) = fuyix) + fdv{^)i sea density S{x) and 
the gluon G{x) for an quark cluster (A^ = 3, 6), which we used in our analysis presented 
here, are given by 



2A^-2 



- 1 

'^'^^^'^^ = 2{AN - 3) + ~ "^^^ 
= xG^{x) = {CN + 1) (1 - xY- , 



with = 9, ag = 11, C3 = 7, and cq = 10 11 . Here B is the usual Euler function and 



Sn{x) represents the sum of sea quark densities over all flavours. With a further assumption 
of f-s, ,n{x) = 2 fd,N{x), the u distribution for the A^-quark cluster is given by 

fu,N{x) = I Sn{x). 
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III. THEORETICAL FORMALISM 



The problem of heavy quark- antiquark pair (QQ) production has been studied exten- 
sively both experimentally and theoretically. In the pQCD approach the quarkonium pro- 
duction cross-section can be calculated by convoluting the hard scattering subprocess cross- 
section (at a give order in a^) with the appropriate initial parton densities and suitably 
chosen hadronization functions. Two popular hadronisation models which describe conver- 
sion of the QQ pair into quarkonia are i) the semilocal-duality (SL) model and ii) the colour 
singlet (CS) model. Since we employ both of them in our work, we will discuss them briefly 
below. As we shall see, the differences in their details result in different kinematics for the 
same process at formally the same order of perturbation theory. It is therefore interest- 
ing to compare their predictions with the E772 data to check the robustness of the pQCD 
approach; any differences are likely to provide a clue on the hadronisation of the {QQ)-peiir. 



A. Semilocal Duality Model 

In this model one first computes the basic QQ production cross-section. The quarkonium 
cross-section is then obtained by simply restricting the invariant mass of the QQ pair between 
2Mq and 2MhQ (where Mq is the heavy quark mass and Mhg is the mass of lowest lying 
g-flavoured meson) and by multiplying the cross-section by a constant. The underlying 
assumption is that irrespective of its invariant mass the {QQ)-peiiT is equally likely to turn 
in to a given quarkonium provided it is below the threshold of the heavy flavoured mesons. 
The constant may depend on the colliding energy. However, we will need to consider only 
ratios of cross-sections for the nucleon and nuclear target at a fixed energy, and the constant 
drops out of such ratios. The QQ cross-section itself is given by, 

fpi/p{Xl)fp2/A{x2) + fp^/p{Xi)fp^/A{x2) 

Pl,P2 " s " sxi 

ap,p,{s,Ml) (8) 



aipA ^ QQX) 



4M^ dx2 
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Here fp/h{x) is the parton density for parton p in hadron h, x is the momentum fraction 
of h carried by p and a is the appropriate subprocess cross-section integrated over all the 
subprocess variables. The square of the total centre of mass (cm) energy of the partonic 
system is given by 

S = SX1X2 

where s is the square of the cm energy of the p-A system. 

For QQ production both the O(a^) Q and O(a^) 
The subprocesses which contribute at 0{a'^) are 

99 ^ QQ, qq QQ (9) 

and the ones contributing at O(a^) are 

99 QQ9, qq QQ9, q9 QQq, q9 QQq (10) 

The quarkonia produced from O(a^) subprocess can have only small pt, being of the order of 
the intrinsic transverse momentum of the incoming partons. On the other hand, the 0(0;^) 
processes will yield large (~ O(GeV)) pt for the quarkonium with the light parton in the 
final state constituting a jet . Thus the J/ip{T) produced at large-pT in the SL model come 
from 2^3 subprocesses. 

The matrix elements squared, appropriately averaged, for these processes have been com- 
puted. There are clearly two classes of diagrams — those involving one gluon, and those with 
three. We shall denote the squared matrix elements by A and B respectively. They depend 
on the momenta of all the five particles. Instead of reproducing the lengthy expressions for 
them, we refer the reader to the original work ||32|. For computing pj^-distributions of the 
quarkonia we will need cross-sections more differential than the expression in eq. (^. 

The phase-space for these processes is 6- dimensional. We consider the kinematics in the 
EHLQ conventions, modified to account for massive Q and Q. Defining the parton 
momenta in their center of mass, one sees that the three final state particles lie in a plane. 

9 



21] expressions for a are available. 



We choose this to be the xy-plane, and orient our axes by choosing the jet to be in the 
positive x-direction. The 4-momenta of the Q, Q and the jet are given, respectively, by 



V s 

P3 = ^ xs{l, (3s COS 635, (33 sin 635,0), 



V S 

Pi = ^ a;4(l,/34 COS 6^45, /34 sin 6*45,0), (11) 



V S 

P5 = ^ a;5(l, 1,0,0), 



where /3j = y 1 — AMg/xfs and the angles are given by 

n {(^AX.y - {(33X3)' - Xj {(33X3)' - {(3,x,r - xi 

co«^35 = ^^5-— and COS ^45 = ^^T— • (12) 

IP3X3X5 /P4X4X5 

Furthermore, energy-momentum conservation implies X3 + X4 + X5 = 2, (0 < Xj < 1), leaving 
only two independent variables in the equations above. In addition to xi and X2 of eq. (||), 
two more variables are needed to specify the momenta of the incoming partons in this frame. 
Choosing these to be the Euler angles {6, 0), we get 



V S 

Pi = - — (1, — sin 6* COS 0, — sin6' sin 0, — cos( 
2 



V S 

P2 = (1, sin 6* COS 0, sin ^ sin 0, COS 6*). (13) 

We can eliminate 6 in favour of the transverse momentum of the pair (which is equal to that 
of the jet) through the relation 



Pt = cos^ 6^ + sin 6' sin 0. (14) 

Similiarly we can write the Feynman scaling variable xp for J/ip as 

xp = - [(xi + ^2)^5 sin6' COS0 + (x2 — Xi)(a;3 + X4)] . (15) 

The E772 data has cuts on this xp which we include in our computations. 

Using the above kinematic relations the fully differential cross section for J/ ■i/'-production 
computed to order can be written down as below: 
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dxidx2dx'idx4^d(j)dpT 1QtisxiX2x\ cos 9 cos^ 

\ H {fq/pi.Xl)fq/A{.X2) + fq/p{Xi)f^lA{X2)] A{p^,Pi, -pi, -p2, P5) 

~^II{/9/p(^i)/9M(^2) + /g/p(a;i)/g/A(a;2)}A(p3,p4,P5,-Pi,-P2) (16) 
'^Yl {f9/pi^i)fq/A{x2) + /g/p(a;i)/3M(a;2)} ^(P3,P4, -Pi,P5, -P2) 

+ ^f9/p{xi)fg/A{x2)B{p4,P3,P5, "Pi, -^2) 

This is the cross section formula which together with appropriate structure functions, cho- 
sen from section II, we integrated over the kinematic region corresponding to the E772 



experimental cuts [lT7| , p!8| to compute the pt distributions for the quarkonia J/ip and T. 



B. Colour Singlet Model 



The colour singlet model was first developed for photoproduction of quarkonia ||3J] and 
later generalised to the hadronic production In this model one projects out the state 
with appropriate spin, parity and colour assignments from the full QQ production amplitude 
to match the quantum numbers of the resonance under consideration. The projection is 
done at the level of hard scattering amplitude itself and this yields a multiplicative factor 
related to the quarkonium wave function at the origin in co-ordinate space. The effect of the 
hadronisation of the QQ pair into the quarkonium is thus contained in this factor. For an 
S-wave resonance this multiplicative factor is the wave function, i?(0)2, at the origin whereas 
for a P-wave resonance it is the derivative of the wave function at the origin, R[{0)^. R{0) 
is related to the measured, leptonic ^5*1 width by 

r«-('5i) = — (17) 

where a is the electromagnetic coupling, eg is the quark charge in units of proton charge 
and M is the mass of the quarkonium. -R'i(O) can be related to the total hadronic width of 
the resonance by assuming it to be approximately the same as its gluonic width given by, 
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r„('Po) ^ (18) 

with as being the running strong couphng. 

The model is known to give a good description of the kinematical distributions in lep- 
toproduction and hadropro duct ion l2^j3^ of J/ip. However, there is a considerable 
uncertainty in the overall normalisation. Even after using the QCD corrected version of 
eq. (p!?!) the data required a /T-factor of 2.4. This large i^-factor is perhaps due to 
the nonrelativistic treatment of the quarkonium J / if) in arriving at the hadronisation factor 
i?(0)^. Once again, for ratios of cross-sections at the same colliding energy, which we will 
consider here, the precise value of ii'-factor plays no role; we assume it to have no nuclear 
dependence. 

The quarkonium production cross-section for the "^^^^Lj quarkonium state in this model 
is given by 

oo I 1 ^max max , max 

fpi/p{xi)fp2/A{x2) + fp2/p{xi)fpr/A{x2) 

xda/di{p,+p2^ ^^+^Lj+ps) (19) 

Here xi,X2 again denote the momentum fractions of the proton and nuclear target re- 
spectively carried by the partons Pi,P2,fpi/h are the parton density distribution functions 
and da/dt denotes the differential subprocess cross-section. Here yi,y2 denote the rapidi- 
ties of the quarkonium and the jet respectively and pt is the transverse momentum of the 
quarkonium. The variables Xi , X2 are given in terms of yi , y2 by 

1 



1 

Xo = - 

2 



xxe + xtc 



(20) 



In the above equation xt = xt = yx^^ + Ar with r = M'^/s. The total allowed range 
of integration over yi,y2, for a given value of px is given by 
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\yi\ < cosh ^ ( — 

V XT 

-in ^^-^TeM-y,)^^ ^^^^^^ ^-xreMvi)^ (2,) 
s, t, u are the Mandelstam variables of the subprocess given by 



s = X1X2S; t = M — Xiy s{pt + M'^)e u = —Xipry/se 

At the lowest order in the strong coupling, O(a^), there is no parton in the final state; 
one has only the gluon fusion process 

gg ^ 'So, 'Po,2 (22) 

The quarkonium so produced has pt — 0. At O(a^), the gq{gq) and qq scatterings give rise 
to ^5*0, ^Pj resonances (x states) via 

gq(q)^'So, 'Pj + qiq); 

qq^'So,'Pj + g (23) 
while the ^Si resonance can be produced directly as well via the gg subprocess, 

gg^ 'So, 'Si, 'Pj + g (24) 
The ^Pj states can decay into "^Si states via the radiative decay 

'Pj^ 3^1+7 , (25) 

thus giving to rise to indirect contribution to the production of ^Si in addition to that of 
eq. ( p^ . The expressions for the differential cross-section for the various subprocesses are 
available |2^J3^ and are not reproduced here. 



Note that according to this model ^Pi production is not possible at O(a^) and hence 
the total ^Pi rates are expected to be small compared to that of ^Po,'^P2. However this 
seems to be belied by the data |^ on hadroproduction of ^Pj charmonium states. Since 



soft gluons, ignored as a part of the model ansatz, are likely to be more relevant in this case, 
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this possibly points to problems with the CS model for the pt — quarkonium production. 
However, we use it here only for large-p^ (~ O (GeV)) quarkonium production where the 
soft gluons most likely do not play such a significant role. 

We compute the ^Si J/ip{T) production cross-section including both the direct con- 
tribution of eq. ( P^ and the contribution from the decay of ^Pj states where ^Pj state 
cross-sections are computed using eqs. (^) and (p4|) . Since the T produced in process in 
eq. (p5|) carries very little energy, we assume the kinematical variables for the ^Pj and ^5*1 
resonance to have essentially the same values. The da/dpT is obtained using eq. ([19|). The 



E772 ||T7| , p!8| cuts on the Feynman xp of the resonance were implemented by restricting 
the ?/i-integration between y™™ and y™^^, which are related to the allowed range of xp, 
Xp <xf < xj., by 

/ 



yr 



max 
ill 



in 



in 



XT 



\ 



-\ 2 



XT 



+ 1 



(26) 



Tables ^ and ^ give the details of the set of the resonance parameters which we used for 
our computations. The masses of all the resonances are taken from the latest compilation 
of particle properties [§. Note here that the value of the i?(0)^ for the S-wave resonances ijj 
and ip' are determined from the updated measurement of the letonic decay width ||^ while 
the R{Oy for the T system and -R'^(O)^ values for both the ip and the T system have been 



taken from ref. 23 



C. Comparison of the Hadronisation Models 

From the above discussion, it should be clear that the two hadronisation models differ 
quite significantly in the time scale for quarkonium formation. In the CS model it is the 
same as the perturbative time scale (~ 1/2Mq ~ 0.07 fm for Q = c) whereas in the SLD 
case it is a typical hadronic scale (~ 1 fm). As a result, the only effect that a medium like 
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the Quark-Gluon-Plasma (QGP) can have in the first case will be on the propagation of 
the quarkonium, whereas in the SLD picture the formation process itself can be affected by 
the QGP environment. 

The kinematics of the two models is also quite different. In the CS case basic hard 
subprocess which gets convoluted with the parton densities (as in eq. (]TP|)) is a 2 — > 1 
subprocess for 0{a'j.) case and 2 — > 2 subprocess for O(a^) case of large-pT quarkonium 
production. In the SLD case, these are respectively a 2 2 subprocess and a 2 — 3 
subprocess. Hence, in principle the momentum division between the large-pT quarkonium 
and the light parton jet can be quite different in the two cases and it is worthwhile to find 
out whether perhaps the px distribution da/dpx for quarkonium could help discriminate 
between these two models of hadronisation. 



D. Gluon dominance 

In Ref. P it was argued that the E772-data are described by 0(al) partonic cross- 
section and further that this production cross-section of eq. (^, is dominated by the 
99 ~^ QQ contribution of eq. Then the ratio of experimentally measured cross-sections 
for different targets, at a given value of Xi;' and r directly yields the ratio of gluon densities 
for the two targets at an X2 given by 

1 



xi,2 = -(xF±^Jxl + r^j . (27) 
It is worth noting that only for the case of 2 ^ 1 kinematics the xp of the quarkonium 



is related to Xi and X2 through the oft-used simple relation of eq. (pT]) . As can be seen 
from eq. ([T5|), this is no longer true for quarkonium production at large-pr via the 2 3 
subprocess of SLD or also the 2 — * 2 subprocess in CS-model. Considering the rather large 
PT-range of the E772-data (< 2.25 GeV for the J/ip and < 4 GeV for the T), it seems to 
be an oversimplification to employ either the eq. ( pTj ) or indeed the O(a^) subprocess itself 
(unless one tolerates rather large values of intrinsic for the partons) . As mentioned 
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above, since the CS-model is perhaps unrehable for ~ 0, even if eq. ( pTl) were to be vahd 
one still further needs the assumption of gluon dominance if the SLD-model of hadronisation 
is employed. 

In support of their argument of gluon dominance Ref. considered the ratio of parton 
luminosities, 

^ ^E31= ^9=M,^ {fq/h{Xi)fq/h{X2) + /g/fe (j^l) /g/fe (3^2) ) 

fg/h{Xl)fg/h{x2) 



for a particular parametrisation of the parton densities in the proton viz. DFLM |^ and 
showed it to be rather small in the Xi^-range of interest (using eq. (|27|) ). It turns out, how- 
ever, that this demonstration of gluon-dominance strongly depends on the choice of parton 
densities used. We plot in Fig.|| the ratio in eq. (pH]) as a function of xp for y/r = 0.0775 
corresponding to the nucleon-nucleon cm. energy of the E772 experiment and Mqq = Mjj^. 
Xi and X2 for a given value oi xp are obtained from eq. (^). The different curves of Fig.|^ 
correspond to various popular choices of the parton densities: the older parametrisations 



DFLM DOl, D02 g and GHR gT[, and the newer MTl, MT2, MT3 |2| and ON [|3 



For the DFLM parametrisation which the authors of Ref. ||^ use the approximation of gluon 
dominance is indeed good for low values of xp {xp ^ 0.3) . However, E772 data goes upto 
Xi? ~ 0.65. Moreover, we also see from Fig.^that gluon dominance crucially depends on the 
choice of parton density parametrisation in a proton. For the MT parametrisation, e.g., the 
qq contribution is ~ 60% of the gg contribution at the largest Xp value considered {xp = 0.5). 
As Fig.^ shows the situation becomes much worse for T with y/r = My/yfs = 0.2365. Thus 
it is clear that the extraction of Pg{x) by using gluon dominance and then simply taking the 
ratio of experimentally measured J/ip-cioss sections for the nuclear and the nucleon target 
has an intrinsic uncertainty of ~ 10-15 %. In view of the fact that the expected deviation 
of pg from unity, viz. the EMC effect, is also of the same order of magnitude, we feel that 
this is not a very effective or precise determination of pg. 

The models for EMC effect always give a parametrisation for nuclear parton densities 
for a specific parametrisation of parton densities in proton. In all the models we use, the 
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choice of reference parton densities in the proton is different from all the above mentioned 
parametrisations. We have checked that our observation about lack of the gluon dominance 
holds for these parametrisations as well. 

In view of the above discussion, attempts to extract the nuclear gluon density from the 
high statistics E772 data appear to have both conceptual problems and sizeable theoretical 
uncertainties. An alternative, albeit a less ambitious, approach may be to check the consis- 
tency between the predicted ratios of differential cross- sect ions da/dpx for various models 
of nuclear parton densities and the high statistics E772 data. One may thus hope to pin 
them down or even expose their deficiencies using the high quality and more differential 
data. Since the small-p^ region is plagued by the issues of intrinsic transverse momentum of 
the inital partons or resummations of higher order diagrams, we choose to investigate only 
the large-pT E772 data and compute only the O(a^) contributions to the cross-sections. 
This compels us to ignore the data on py-integrated xp distributions or their ratios since 
they should receive substantial contribution from the small-pr region. In order to test the 
Xi?- dependence of our pQCD calculations, it would be desirable to have the data integrated 
in different, at least two, ranges of pt- 



E. Results 

The E772 experiment has provided data for the ratio 

dpT I dpT 

with an xp-cut of 0.15 < xp < 0.65 on the J/ip's, while for the T-production cross sections, 

they chose to present only a{pT), where 

dajpA TX) ^ ^^^(p^) da{pp TX) 
dpT dpT 

with a corresponding xp-cut for T of —0.2 < xp < 0.6. The nuclei used were carbon, 

calcium, iron and tungsten. Incorporating these XiT'-cuts, using eqs. (^) and (^), we com- 



puted each of the individual pr-distribution in eqs. ( P9| - |30|) for both the SLD-model and the 
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CS-model and the three models of nuclear structure functions discussed in sect. II. Fig. ^ 
exhibits our results for the SLD-model for all the four nuclei along with the corresponding 
data from the E772 collaboration. The errors for theoretical predictions are purely statis- 
tical, arising from the Monte Carlo integration of the differential cross sections. One sees 
that for the lighter nuclei both the two-component models, namely, the gas model and the 
six-quark cluster model, describe the data well, especially if one takes into account a possible 
systematic error of a few per cent due to variations in input parameters such as Mc, Aqcd 
etc. For the tungsten nucleus, however, none of the models seems to to be in agreement with 
data. One could compare the CS-model results with the E772-data in a similar manner 
as well. Instead we choose to compare the results of the two hadronisation models directly 
since the theoretical results have a better precision. Fig. |^ shows the ratio of R^/''^{pt) for 
the SLD-model and the CS-model for each model of the nuclear parton densities, nucleus 
and pt- The ratio of the ratios seems to lie in almost all the cases within ~ 5% of unity. 
Considering the different kinematics of the two models and also the different physics of the 
quarkonium formation the agreement is really remarkable and it shows the robustness of the 
pQCD predictions. Note that even in the case of Tungsten both the models agree with each 
other rather well and thus have essentially the same discrepancy with the E772 data. Of 
course, the discrepancy in the case of tungsten does expose the inadequacy of all the three 
models of the EMC effect and the corresponding parametrisation of the nuclear parton den- 
sities considered here but the generalagreement in other cases, on the other hand suggests 
that the structure function effects indeed do describe the bulk of the p^-data. One may 
note here that the lowest pt value at which we performed these computations is somewhat 
low, being 0.79 GeV. Presumably, the anticipated large QCD-corrections at such low px 
affect all the cross sections similarly and thus cancel out in the ratio, resulting in a good 
description of the data. A remark on the gluon dominance may be in order as well. Extract- 
ing the contribution for individual parton subprocesses, we typically found the dominant 
contribution to be from gg and qg processes which were in the range of 75-80% and 25-20% 
respectively. Thus, the quark contributions are sizeable and an independent determination 
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of the gluon density from these data are not possible without making an ansatz for the quark 
distributions. 

Fig. ^ shows a comparison of our calculations for the SLD-model for T-production with 
the E772-data. At each pt in the range of 1-4 GeV, the calculations for each individual 
da/dpT were done as for J/^/^-pro duct ion by incorporating the experimental Xi^-cut. The 
results for differential cross sections were then fitted as a power law in A to obtain a{pT)- 
One sees a similar general agreement for the gas model and the six-quark cluster model as 
for J/'?/'-production at moderate values of pt- At the largest pt, however, the E772-data 
rise too sharply compared to any model and could possibly indicate that these models tuned 
to earlier large x-data have to be better tuned to perform well in the small x-region. Fig. 
1^ shows that this disagreement at large pt, as well as the agreement at lower p^, are once 
again features which do not depend on the hadronisation model. The figure shows the ratio 
of the results obtained for the SLD-model and the CS-model for the experimental cuts of 
E772 as a function of pt and they are within ~ 5% of each other. Thus within the framework 
of pQCD these data too are explicable in terms of changes of nuclear structure functions. 
As one can expect, the contribution of the qg subprocess in the case of T-production is 
even larger, being typically ~40% in comparison to that of the gg subprocess which almost 
accounted for the rest. 



IV. DIMUON PRODUCTION 

The J/ip{T) is detected most efficiently via its decay into a fi^fi^ pair. Hence any critical 
evaluation of J/ ip suppression as a signal of QGP formation also needs a good understand- 
ing of the Pt- and A-dependence of this continuum fi'^fi~ background. An experimental 
measurement of the A-dependence of the dimuon pair production can also provide an inde- 
pendent probe of the nuclear parton densities and an evidence in favour of the universality 
of the EMC-effect, i.e., its process independence. In fact, one of the earliest theoretical 
attempts to understand the dimuon data had postulated |Q A-dependent sea density even 
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before the EMC effect was experimentally discovered. The production of massive fi'^fi~ 
pairs (DY) in hadronic collisions is now well understood in the framework of pQCD ||22|| . 
The dimuon production at small pt and large is essentially well described in terms of qq 
annihilation process in spite of the large higher order corrections. For large pt of the fi~^fi~ 
pair the production cross-section is given by the 0{as) subprocesses involving gluons viz., 

q + 9-^l* + q^ +g; q + q-^-f*+g^ + 9 (31) 



The high statistics E772-experiment [|T^] has provided data for nuclear-dependence of proton- 
induced pair production over a wide range of xp and px values. The data on the ratio of 
the integrated dimuon yield for different nuclei were compared with theoretical predictions, 
obtained by using the qq annihilation process, for various models of the EMC effect. It 
seemed to rule out the 6-quark cluster model |2^. However, a later comparison |Tl 



with an improved version of the model, described in sect. II. C, showed that this model too 
can be consistent with the information on the ratio of the integrated dimuon yields. The 
E772 experiment [|Ty| has also presented pt distributions (integrated over xp and M^+^_) 
and Xp distributions (integrated over pp and M^+^_) for the dimuon pairs. A comparison of 
ratios of these differential distributions for different targets with the predictions of various 
models of EMC effect can discriminate between them more effectively. Since the 0{as) 
pQCD calculation is valid only at large-p^, we restrict ourselves to the pr distributions. 
The xp distributions are integrated over the complete range of pp and hence dominated by 
Pt ^ data, which once again forced us to ignore them in this leading order pQCD analysis. 

The kinematics of the DY pair production at large pp in the 0{as) subprocesses 

of eq. (^) is very similar to the kinematics of the J /ip{T) production in the colour-singlet 
model, discussed in sect. III.B. Its differential cross-section is given by 

da /-s/™ rvT'''' r da{qg ^ fi^fi^q) 



ry'i ry2 r 

/ dyi / dy2 xiX2\ Pc 



where 



20 



Pqg = [fq/pi^l)fg/A{X2) + fg/p{Xl)fq/A{X2) 

q,Q 

Pqq = [fq/pi.^l)fq/A{.y2) + fq/p{.Xi)fq/A{,X2)\ , (33) 



with 



da 



{qg fi^fi q) 



di " 9M2 -rs^u 



Here yi,y2 are the rapidities at the fi~^fi~ pair and the associated jet respectively. The 
Mandestam variables s, i and u, the relation of the momentum fractions Xi and X2 in terms 
of yi,y2, the integration limits and their relation with the experimental Xi^-cut [xp > 0) 
are precisely the same as those given in sect. lll.B for the colour singlet model with 
replaced in all the formulae by M^+^_ . We will therefore not repeat them here. 
Experimental information [|19| is available for the ratio 

i?^^ = ^ (M - /iV-x) / ^ (PP - /^V"^) (35) 

where da^^ /dpr is the differential DY cross section integrated over the continuum region 
(avoiding the resonances) 4 < M^+^- < 9 GeV and M^+^- > 11 GeV, with xp > 0. 

Integrating eq. (^) over the above experimental cut of Xi;' > and 4 < M^+^- < 9 as 
well as M^+^- > 11 GeV, we compute the ratio R^"^ of eq. ( P^D for all the different nuclear 
targets used in the experiment for each of the three sets of nuclear and nucleon parton 
densities described before. M'^+^- was used as the scale for as in eqs. (|34D. 

Fig. 1^ exhibits the results of our computation for the four different nuclei with the 
corresponding data. Again we see similar to the case of resonance production that the 
general trends of the data are well described by the model predicltions for the gas model 
and the 6-quark cluster model. 
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V. CONCLUSION 



In conclusion, we have shown in this paper that the high statistics data E772 on the nu- 
clear dependence of the production of quarkonia (J/^E' and T) and dimuon pairs at large Pt 
can be entirely explained in terms of the same nuclear structure functions, in the frameowrk 
of pQCD. All our theoretical calculations contained no arbitrary free parameters; only ex- 
isting models of the EMC-effect with their already fixed values of parameters were used. We 
employed two popular models of hadronization of the QQ-paii into the quarkonia. In spite 
of their big kinematical differences, we found both to yield predictions which were within 
~ 5% of each other. This shows the robustness of the pQCD approach and underlines the 
importance of the nuclear structure function effects in understanding the behaviour of these 
data. Similar conclusions |^,|^ about the independence of hadronization mechanism and the 
universality of the nuclear structure functions in various hard scattering processes have 
been obtained before but the accuracy of the present data makes them now much stronger. 
Recently it was argued that quantum mechanical coherence and interference effects de- 
stroy factorisation in quarkonia production and hence prevent the possibility of using the 
same nuclear structure functions for different final states. However, the consistency of both 
the quarkonia and the dimuon data with our calculations, points towards the correctness of 
ideas of universality of the structure functions, at least in the kinematic region probed in 
our analysis. 

Our analysis also indicates that the accuracy of the data has now reached a stage so as 
to distinguish between various models of the EMC effect and the nuclear structure function 
parametrisations therein. Indeed, the inability of all models to make a better prediction for 
the tungsten nucleus than shown in Fig. 4, and the disagreement at the largest pt value in 
Fig. 6, are hints for inventing better parametrizations of the nuclear dependence of the quark 
(and gluon) distributions. We also argued that the twin assumptions of gluon dominance 
and adequacy of the lowest order partonic cross section are unreliable due to the large xp 
and Pt ranges of the E772-data. Extraction of the nuclear gluon density using them is likely 
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to be dominated by uncertainties as large as the gluonic EMC-effect itself. 

We have examined here data at a fixed value of pt but which have been integrated over 
the entire xp region corresponding to the acceptance of the experiment. The integrated 
data are dominated by the data at small xp values ( or not-so-small X2 values). Hence our 
non-inclusion of any shadowing effects for the nuclear parton densities can be justified. If 
one wants to critically use these data to study the shadowing effects in the nuclear parton 
densities, then it would be necessary to look at the nuclear effects in the pr integrated data 
at large values oi xp (which will probe small values of X2). The data on xp distributions 
available currently is integrated over the entire range of pr whereas our pQCD analysis is 
valid only for pr > i GeV. If information about the xp distributions integrated over only 
the large pt region becomes available, it will help unravel the issue of nulear dependence of 
the J/^*, T and the dimuon production further. 
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TABLES 

TABLE I. Model parameters for nuclei used in E772 experiment for gas model(T,a;) and six 



quark cluster model(€). 


A 


T (MeV) 








e 


12 


54 




0.069 




0.112 


40 


47 




0.057 




0.170 


56 


45 




0.117 




0.186 


184 


42 




0.132 




0.230 




TABLE II. Resonance parameters used foi 


' 'il'. ih 


and T.T'. 




Resonance R 


M(MeV) 






BF{R- 


-> ^Si+ neutrals) 




3096 


0.542 






1.00 




3686 


0.307 






0.55 


T 


9460 


4.54 






1.00 


T' 


10020 


2.54 






0.19 


TABLE III. Resonance parameters used for x"^ and states. 


Resonance R 


M(MeV) 


i?i(0)VM2(Gey)3 




BF{R ^3 Si+ neutrals ) 


Xo 


3415 


9.1 X 10-3 






6.6 X 10-3 


xf 


3510 


9.1 X 10-3 






0.273 


X2 


3555 


9.1 X 10-3 






0.135 




9860 


1.5 X 10-2 






0.040 


X? 


9890 


1.5 X 10-2 






0.290 


X^2 


9915 


1.5 X 10-2 






0.220 
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FIGURES 

FIG. 1. Predictions for the ratio of the nuclear and nucleonic gluon density, pg, of the three 
models of the EMC effect described in the text. 

FIG. 2. The ratio Rqg = of Eq. (28) corresponding to J jil) production at the FNAL ener- 
gies for MTl(solid), MT2 (long-dashed), MT3(medium-dashed), ON(short-dashed), GGR(dotted), 
GHR(dot-dashed) , DO 1 (dot-dot-dashed) and D02 (inverted triangle) parametrisation of the par- 
ton densities. 

FIG. 3. Same as in Fig. Q but for Upsilon production at the FNAL energies. 

FIG. 4. E772 data on the ratio R'^l'^ of Eq. 29 compared with the predictions for the gas 
model(squares), six-quark cluster model(circles) and the rescaling model(open circles) of the EMC 
effect, obtained using the SLD for hadronisation. 

FIG. 5. Ratio of the predictions for 'R^l'^ for the SLD and CS model of hadronisation, for the 
different nuclear parton densities. Notation is same as in Fig. |^. 

FIG. 6. E772 data on a(pT) of Eq. 30 compared with predictions of the three different models 
of the EMC effect mentioned in the text, using the SLD model of hadronisation. Notation is same 
as in Fig. 

FIG. 7. Ratio of ol{pt) predicted in the SLD and CS models of hadronisation for the three 
models of the EMC effect. Notation is the same as in Fig. ^. 

FIG. 8. E772 data on the ratio of eq.^ compared with predictions of the three models of 
the EMC effect. Notation is same as in Fig. |^. 
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